Efimov physics in bosonic atom-trimer scattering 
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Bosonic atom-trimer scattering is studied in the unitary limit using momentum-space equations 
for four-particle transition operators. The impact of the Efimov effect on the atom-trimer scattering 
observables is explored and a number of universal relations is established. Positions and widths of 
tetramer resonances are determined. Trimer relaxation rate constant is calculated. 
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PACS numbers: 31.15.ac, 34.50.-s, 34.50.Cx 



Few-particle systems with resonant interactions, char- 
acterized by the two-particle scattering length a being 
much larger than the range of the interaction, were 
predicted to have a number of universal (interaction- 
independent) properties and correlations between observ- 
ables [1-6]. Such a behavior was recently confirmed in 
cold atom physics experiments 0, H|, but can be seen 
qualitatively also in few-nucleon systems [l|, 0. One 
of the best-known examples is the existence of the in- 
finite number of weakly bound three-boson states (Efi- 
mov trimers) in the unitary limit a = oo [![. In that 
limit, depending on the available energy, an infinite num- 
ber of atom-trimer channels may be present in the four- 
boson system. Due to complexity of the multichannel 
four-particle scattering problem the universal properties 
of such a system, i.e., the atom-trimer continuum, are not 
known yet; we therefore aim to study them for the first 
time in the present work. Furthermore, the existence of 
a pair of four-boson states (tetramers) associated with 
each Efimov trimer was predicted [H, However, only 
the two tetramers associated with the trimer ground state 
are true bound states that have been studied in all the 
details using standard bound state techniques 0-H, Q- 
Higher tetramers are resonances above the atom-trimer 
threshold and for this reason their properties are far less 
known; they will be determined in the present work using 
proper scattering calculations. 

Our description of the four-boson scattering is based 
on the Alt, Grassberger and Sandhas (AGS) equations 
10] for the transition operators; they are equivalent to 
the Faddeev-Yakubovsky equations [ll[ for the wave- 
function components. Symmetrized form of AGS equa- 
tions [l2| is appropriate for the system of four identical 
bosons, 

Un = Pzi(G$tGo) 1 + P^UiGotGoUn + U^GotGoUzx, 

(la) 

U 21 = (1 + P^iGotGo)- 1 + (1 + P^UiGotGoUn, 

(lb) 

where Up a are the four-particle transition operators, Go 
is the four free particle Green's function, and P34 is the 
permutation operator of particles 3 and 4 that ensures 
correct permutation symmetry of the system. The dy- 
namic input is the two-boson potential v from which the 
two-particle transition-matrix t and the AGS transition 



operators U a are derived, with a — 1 and 2 correspond- 
ing to the 1 + 3 and 2 + 2 subsystems, respectively. As 
explained in Ref. [T2I ] . the atom-trimer scattering ampli- 
tudes are given by the on-shell matrix elements of Lin 
calculated between the Faddeev amplitudes of the corre- 
sponding initial and final states. 

We solve AGS equations using momentum-space 
partial-wave representation [12| where they are a sys- 
tem of coupled three-variable integral equations that af- 
ter the discretization of momentum variables becomes 
a very large system of linear algebraic equations. In 
the case of the four-nucleon scattering those equations 
have been successfully solved with realistic nuclear and 
Coulomb interactions [UlEi]]- In the present calculations 
we take the numerical techniques for the treatment of 
four-particle permutations and trimer bound-state poles 
from Refs. [12Hl4j ]. However, an important difference as 
compared to the four-nucleon system is the presence of 
sharp resonances in the four-boson system, manifesting 
themselves as poles of the AGS operators (Q}. Since the 
convergence of the multiple scattering series in the vicin- 
ity of the pole is very slow, in this work we solve sys- 
tems of linear equations by the direct matrix inversion 
instead of the iterative double Pade summation method 
[l2l ]. Since we are interested in the universal properties 
that must be independent of the interaction details, we 
use rank 1 separable two-boson potentials v — \g)X(g\ 
acting in S'-wave only and thereby reducing Eqs. ^ to 
a small system of two-variable integral equations. Al- 
though the two-boson interaction is limited to S'-wave, 
i.e., l x = in the notation of Ref. [12T |. higher angular 
momentum states with l Vl l z < 2 have to be taken into 
account for the relative motion in 1 + 2, 1 + 3, and 2 + 2 
subsystems to achieve the convergence. 

Our standard potential has simple gaussian form fac- 
tor, i.e., its momentum-dependence is (k\g) = e~( fc / A ) . 
The strength A is chosen to reproduce infinite two-boson 
scattering length. In that limit all observables scale with 
A; e.g., the binding energy of the n-th excited trimer 
b n ~ A 2 . It therefore makes no sense to specify partic- 
ular value of A as well as boson mass. Instead, we will 
use dimensionless ratios. As the length scale associated 
with the n-th excited trimer we will use l n = h/^/2^ib n 
where fi± is the reduced atom-trimer mass. 

We do not include explicit three-body force, however, 
many-body forces are simulated by a different off-shell 
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515.035 


22.695 


4.06 x 10~ 8 


5 


22.6 


-1.09 


3.22 


-0.017 


1 


2126.36 


11.349 


2.14 x 10 -3 


1 


11.7 


-0.21 


2.93 


-0.012 


2 


518.570 


21.528 


8.87 x 10" 5 


2 


22.5 


-1.52 


3.22 


-0.024 


3 


515.042 


22.655 


3.90 x 10 -6 


3 


22.7 


-1.09 


3.22 


-0.017 


4 


515.035 


22.694 


1.72 x 10 -7 


4 


22.6 


-1.09 


3.22 


-0.017 



TABLE I. Trimer properties obtained with one-term (top) 
and two-term form factor potentials (bottom). 



TABLE II. Atom-trimer scattering length and effective range 
obtained with one-term (top) and two-term form factor po- 
tentials (bottom). 



behavior of the two-body potential. To prove that this 
doesn't change universal properties of the four-boson sys- 
tem we also use the potential II with two-term form fac- 
tor (k\g) = [1 + c-2 (fc/A) 2 ]e~( fc / A ) ; large negative value 
of C2 = —9.17 ensures very different off-shell behavior. 
In the configuration space representation the potential II 
has several attractive and repulsive regions much like the 
one of Ref. [l5| and supports a deeply bound trimer that 
is non Efimov state in contrast to excited states and all 
states of the standard potential. This can be seen in Ta- 
ble U which collects ratios for calculated trimer binding 
energies b n : all b n _i/b n are quite close to the character- 
istic Efimov value of 515.035, the exception being &0/&1 
for the potential II that is much larger. However, the Efi- 
mov value is reached with good accuracy only for high 
excited states n > 3. This is not surprising since the 
Efimov condition R n » p ensuring truly universal be- 
havior, where R n is the size of the n-th trimer and p 
the range of the interaction, is only well satisfied for n 
large enough whereas for the ground states Rq < p may 
take place [15j, [16|]. Simultaneously k„_i/k„ converges 
towards 22.694 where n n is the expectation value of the 
n-th excited trimer internal kinetic energy. The ratio 
bn/Kn is a measure for the high-momentum components 
in the trimer wave function; as Table U demonstrates, it 
differs significantly for the two employed potentials. 

In Table ITI1 we present results for the scattering length 
a n and effective range r„ for the atom scattering from 
the n-th excited trimer up to n = 5. For both potentials 
a-n/ln and r n /l n converge towards universal values as n 
increases, i.e., 



^n I In 
Tn f In 



22.6- 1.09?:, 
3.22- 0.017i, 



(2a) 
(2b) 



however, significant potential-dependent deviations are 
observed for n < 2. Including strong repulsive three- 
body force that enforces the Efimov condition R n >> p 
(but with additional numerical complications) probably 
could speedup the convergence with n but even in such a 
case n = would be insufficient since it doesn't account 
for inelasticities. 
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FIG. 1. (Color online) S-, P- and D-wave phase shift for the 
atom scattering from the n-th excited trimer. 



It turns out that the universal limit exists for all scat- 
tering observables. In Figs. Q] - [2] we show all relevant 
phase shifts 5% and S'-wave inelasticity parameter 77^ 
for the atom scattering from the n-th excited trimer as 
functions of the relative kinetic energy E n divided by 
the respective b n ; the elastic S'-matrix is parametrized 
as = Vn e2lS " ■ Again, results with n > 3 are in- 
distinguishable and represent the universal values, but 
deviations are seen for lower n. Elastic scattering is de- 
termined by relative atom-trimer S-, P-, and D-waves. 
In fact, S- wave dominates at lower energies but close to 
E n /b n = 1 the individual contributions of S-, P-, and 
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FIG. 2. (Color online) S-wave inelasticity parameter for the 
atom scattering from the n-th excited trimer. 
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FIG. 3. (Color online) Elastic and inelastic cross sections for 
the atom scattering from the n-th excited trimer. 



D-waves to the elastic cross section are 31, 58, an 11%, 
respectively, while the F-wave with \8^\ < 0.4° yields less 
than 0.1% and therefore is negligible. Situation is differ- 
ent in the inelastic scattering where only S-wave con- 
tributes significantly: only n^ clearly deviates from 1 as 
shown in Fig. [5] while in higher partial waves, due to very 
different size of trimers and the angular momentum bar- 
rier, transitions to other trimers are strongly suppressed 
and n^ are very close to 1, e.g., 1 — r]£ < 10~ 6 . 

The elastic and inelastic cross sections a(n — > n') for 
the atom scattering from the n-th excited trimer leading 
to the lower- lying n'-th trimer are presented in Fig.[3J We 
do not show cross sections for n < n' that can be obtained 
by time reversal as er(n' — > n) = (E n / E n i)a(n — > n'). 



Within the resolution of the plot <r(n —> n')/l„ for 
n, n' > 2 and fixed n — n' are independent of n and 
employed potential and thereby represent the universal 
values. Furthermore, for n > n' large enough the ratios 



n 



a(n — > n' — 1) 



43.7 



(3) 



are energy- independent. These observations allow us to 
extrapolate a(n — > n') results to any n and n' that are 
sufficiently large. 

Perhaps most interesting energy regions, namely those 
containing S*-wave four-boson resonances just slightly be- 
low E n /b n — 1, are not displayed in Fig. [3J There are 
two (k = 1, 2) four-boson resonances associated with the 
iV-th Efimov trimer. Resonances are poles of the AGS 
transition operators in the complex plane; thus, in the 
vicinity of the four-boson resonance 



Up a « U^\E - Er)- 1 +U^> +U^(E - Er), (4) 

where E = E n — b n and E r = —BN,k ~ i^N^k/^ with 
—BN,k being the (N, k)-th resonance position relative 
to the four-body breakup threshold and Tjv^ its width. 
These parameters were determined by fitting the on-shell 
matrix elements of U\\ into Eq. (Q}. Again, for N large 
enough, i.e., N > 3, the results with good accuracy be- 
come independent of potential and N, 



B NA /b N « 4.6108, 
B N . 2 /b N « 1.00228, 



rAr,i/2&Ar 

rAf,2/2&AT 



0.01484, (5a) 
2.38 x 10~ 4 . (5b) 



We note that Bn k but not Tx,k have already been cal- 
culated in Ref. @. While B N<1 /b N w 4.58 of Ref. 5] 
is quite close to our number, the k = 2 resonance with 
BN,2/bN ~ 1-01 was predicted in Ref. [j| to be signifi- 
cantly further from the atom-trimer threshold than our 
result. In contrast to S'-wave, there is no four-boson res- 
onances in higher angular momentum states as can be 
seen in the bottom panel of Fig. [TJ 

Unlike in nuclear physics, the direct measurement of 
the atom-trimer cross sections in cold atom physics ex- 
periments is not possible yet. Instead, one may be able 
to create an ultracold mixture of atoms and excited Efi- 
mov trimers in a trap and observe the trimer relaxation, 
i.e., the inelastic collision of an atom and trimer in the 
n-th excited state leading to the atom and trimer in the 
lower-lying n'-th state. The kinetic energy AK w b n i 
released in this process is shared between the atom and 
trimer with the ratio 3:1. Thus, if b n i/A is larger than 
traping potential, the final-state trimer escapes the trap. 
Assuming that this is the dominating mechanism for the 
trimer loss, the time evolution of the density p n (t) of the 
n-th excited state trimers in the trap is given by 



dp„(t) 
dt 



-(3nPa(t)p n (t), 



(6) 



p a (t) being the atom density and f3 n the relaxation rate 
constant |2j. The alternative way of the trimer loss, 
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FIG. 4. (Color online) Temperature dependence of the trimer 
relaxation rate constant. 

i.e., inelastic trimer-trimer collisions, is suppressed if 
Pn(0) << Pa(0). Under this condition p a (t) ~ Pa(0) 
and Eq. © has a simple solution 

Pn(t) =p„(0) e -^<°)*. (7) 

Thus, in this case the lifetime of the mixture is sim- 
ply given by l//3„p a (0). The relaxation rate constant 
Pn = J2n' fin^n' has contributions /3 n ^, n > = (v n a(n -> 
n')) from transitions to all trimers n! < n, where v n — 
y^2E n / [i\ is the relative atom-trimer velocity and (. . .) 
denotes the thermal average. Thus, the trimer relaxation 
rate constant is determined by the atom-trimer inelastic 
cross sections calculated in the present work. In par- 
ticular, Eq. ([3]) implies that for n and n' large enough 
Ai-s-n' / 'fin^ri -l ~ 43.7 and therefore the relaxation is 
strongly dominated by the n — > n — 1 transition. The 
zero temperature limit of the relaxation rate constant 
can be obtained using the optical theorem as 

/£ = -— Man). (8) 
Mi 

The results at finite temperature T are given in Fig. @] as- 



suming the Boltzmann distribution for the relative atom- 
trimer energy; ks is the Boltzmann constant. Figure 0] 
indicates that the use of T = limit is inappropriate at 
temperatures above ksT/bn > 1CP 4 . 

In summary, we studied bosonic atom-trimer scatter- 
ing in the unitary limit. It is a complicated multichannel 
four-particle scattering problem involving, in the present 
calculations, up to six open channels with the Efimov 
trimer binding energies differing by a factor larger than 
515 5 « 4 x 10 13 . Exact AGS equations were solved in 
momentum-space framework with some important tech- 
nical modifications compared to previous calculations of 
the four-nucleon system. The results for reactions with 
highly excited trimers (at least 2nd excited state) in the 
initial and final channels were found to be independent of 
the used potential and thereby represent universal values 
for atom-trimer scattering length, effective range, phase 
shifts, elastic and inelastic cross sections and four-boson 
resonance parameters. On the other hand, results for 
lower trimers demonstrate that significant quantitative 
deviations from the universal behavior are possible. The 
comparison with the experimental data could not be per- 
formed yet, but the obtained atom-trimer scattering re- 
sults were related to the trimer relaxation rate constant 
that hopefully will be measured in the future experiments 
with ultracold mixtures of atoms and Efimov trimers. 

The developed technique is applicable also to dimer- 
dimer scattering. Our first calculations confirm the 
nontrivial behavior of the dimer-dimer scattering length 
when approaching the unitary limit [4J| but can provide 
results also at finite energies. Furthermore, the exten- 
sion to fermionic systems in the unitary limit may have 
impact not only on the cold atom but also on nuclear 
physics, e.g., by clarifying to what extent four-nucleon 
resonances are universal. 

The author thanks R. Lazauskas and L. Platter for 
discussions and suggestions. 
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